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Abstract 

This paper presents a thorough evaluation of a bistable system versus a matched filter 
in detecting bipolar pulse signals. The detectability of the bistable system can be opti- 
mized by adding noise, i.e. the stochastic resonance (SR) phenomenon. This SR effect 
is also demonstrated by approximate statistical detection theory of the bistable system 
and corresponding numerical simulations. Furthermore, the performance comparison 
results between the bistable system and the matched filter show that (a) the bistable 
system is more robust than the matched filter in detecting signals with disturbed pulse 
rates, and (b) the bistable system approaches the performance of the matched filter 
in detecting unknown arrival times of received signals, with an especially better com- 
putational efficiency. These significant results verify the potential applicability of the 
bistable system in signal detection field. 

1 Introduction 

Stochastic resonance (SR) phenomena involve optimizing a system performance measure 
via adding noise, and have received considerable attention in the past few decades [1-22]. 
The prototype SR model is an overdamped bistable system, which was used to describe the 
earth's climatic change [1-3]. Since then, this generic system is a frequently used model for 
characterizing SR phenomena in diverse scientific fields [4-14]. Moreover, many experimental 
verifications of SR effects have been demonstrated for a bistable characteristic, such as in 
a Schmitt Trigger [15], bistable ring laser [16], paramagnetically driven bistable buckling 
ribbon [17], bistable electron paramagnetic resonance systems [18], bistable superconducting 
quantum interference devices [19], vertical cavity surface emitting lasers [20], etc. In addition, 
new SR-type phenomena were also observed in bistable systems subjected to suprathreshold 
input signals [21,22]. 

Recent research has focused on potential SR applications in signal transmission [20,23-27] , 
estimation [28,29] and detection [30-48]. Signal detection theory offers a powerful tool for 
analyzing linear systems, such as the matched filter with a Gaussian white noise background 
and its impulse response matched to the input signal maximizing the filter output signal- 
to-noise ratio [49]. These nonlinear systems, thereby, are potentially rich in signal detection 
applications, though generally more difficult to theoretically tackle [5,7,34,35]. A variety 
of appealing nonlinear systems or models exhibiting SR effects have been explored as detec- 
tion devices, such as the threshold detector [35,37], electronic circuits [5,40-42], biological 
systems [5-7,10,13], neuron models [45,46,48], bistable systems [30-32], etc. Although 
many researchers stressed that SR can not do wonders in signal detection, in the sense 
that nonlinear systems exploiting SR can never do better than the corresponding linear sys- 
tems [30-32,46,50-52,61] in the presence of Gaussian white noise, some promising results have 
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been obtained for non-Gaussian noise conditions [34-36,43,44], the computational and design 
efficiency of nonlinear systems [32,53], noise-floor limited systems [54,55], psychophysical ex- 
periments [23] or models [56] , and parallel arrays of nonlinear devices via the suprathreshold 
SR effect [57,58,62-64]. 

In line with this, we investigate further the performance of bistable systems versus 
matched filters in detecting bipolar pulse signals. Some new merits of bistable systems will 
be emphasized in the comparison of results. Section |21 develops a theoretical non-stationary 
probability density model of the bistable system. This approximate theory describes the 
SR phenomenon in the bistable system well, and agrees with the numerical simulation re- 
sults. The comparison of the pertinent signal detection statistics, between bistable systems 
and matched filters, demonstrates the fact that the bistable system is a suboptimal detec- 
tor against the matched filter [30-32,57,58]. Beyond this fact, we evaluate in detail the 
detectability of bistable systems, versus matched filters, in detecting pulse signals with dis- 
turbed pulse rates and unknown arrival times in Sec. [21 From the power spectrum analysis 
of bistable system outputs, the different pulse rates of the disturbed signal can be clearly 
discriminated. The matched filter, however, is more like a lowpass filter, wherein too much 
noise at low frequency bands weakens its resolution. In the case of detecting the unknown 
arrival time of received signals, the bistable system, with lower computational overhead, is 
an appropriate substitute for the matched filter. The corresponding numerical simulation re- 
sults also indicate that the bistable system is a more competitive detector than the matched 
filter, with respect to its robustness and computational efficiency. This may be of importance 
to bistable electronic or optical devices utilized in signal detection. We argue that the po- 
tential applicability of the bistable system deserves further study in other signal processing 
problems. 

2 Signal detection statistics of the bistable system ver- 
sus the matched filter 

Here, we consider the detection problem of a bipolar pulse signal s(t) and introduce two 
hypotheses 

H : s{t) = -A, [(n-l)T p , nT p ], 

fTi:a(t) = +A, [(n-l)T p , nT p ], (I) 

where n — 1,2, . . .. The input pulse signal s(t), for (n — f)T p < t < nT p , has the pulse 
duration time T p and pulse amplitudes ±A (see an example of Fig. f^l . The background 
noise r](t) is additive Gaussian white noise with autocorrelation (rj(t)r](0)) — 2D5(t) and 
zero-mean. Here, D denotes the noise intensity. The mixture of s(t) plus rj(t) is then applied 
to a detector, such as the matched filter or the bistable system introduced below. 



2.1 Nonstationary probability density model of the bistable system 

The bistable system under study is described as 

r a ^-=x(t)-^-+s(t)+ V (t), (2) 

with real system parameters r a and Xb [25]. r a is related to the system relaxation time. 
The dynamics of Eq. (J2J) are derived from the symmetrical double- well potential Vo(x) — 
—x 2 /2 + x 4 / (4X 2 ), having the two minima Vb(±Xfc) = — X 2 /4. Parameters r a and Xf, have 
the units of time and signal amplitude respectively, and define natural scales associated to 
the process of Eq. @ . 

In each pulse duration T p , the system of Eq. J5| is subjected to the constant amplitude 
+A or — A, with an additive input Gaussian white noise r/(t). Under these conditions, the sta- 
tistically equivalent description for the corresponding probability density p(x, t) is governed 
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by the Fokker- Planck equation 



dp(x,t) v d T . D d 2 n , 

where V'(x) = —x + x 3 /X 2 =p A and the Fokker-Planck operator is Lpp = -J^V'(x) + ^-j^i 
[65]. Here, p(x,t) obeys the natural boundary condition that it vanishes at large x for any t. 
The steady-state solution of Eq. for a permanent input at +A or —A, is given by 

p(x) = lim p(x,t) = Cexp[-^M], (4) 

where C is the normalization constant [65]. This solution of Eq. Q is for global equilibrium 
in the double-well potential. 

Now, we will seek the nonstationary solution p(x, t) of the Fokker-Planck equation, 
Eq. (0), in the case of an input transition from s(i) = —A to s(t) = +A, or vice versa. 
This calculation is carried out in Appendix A. We show in Appendix A that the transition 
from the stationary density corresponding to s(t) = —A to the stationary density corre- 
sponding to s(t) = +A, or vice versa, is dominated by an exponential temporal relaxation 
with a time constant 1/Aj.. This allows us to deduce a response time T r = T a j\\ for the 
bistable system, which is a measure of the time taken by the system to switch from one 
potential well to the other, when the input changes from ±A to =fA, in the presence of noise. 
In Appendix A, the non-stationary solution p(x,t) for an input transition from s(t) = —A 
to s(t) = +A (or vice versa) is approximated with the two first terms from its asymptotic 
representation of Eq. (|25(l , as 

p(x, t\s(t) = ±A) ~ P {x\s{t) = ±A) + [p(x\s(t) = tA) - P {x\s{t) = ±A)} exp(-i/T r ), (5) 

where p(x\s(t) = ±A) are the steady-state solutions of Eq. (@J. In Eq. J^J, when t = 0, the 
term exp(— t/T r ) = 1 and p(x,t\s(t) = ±A) starts with the initial condition of p(x\s(t) = 
=FA). As t — > +oo, the term exp(— t/T r ) = 0, and p(x,t\s(t) — ±A) tends to the stationary 
condition of p(x\s(t) — ±A). This theoretical model, although approximate, nicely captures 
the double role played by the noise, both in accelerating the switching dynamics between 
wells and in enhancing the fluctuations inside the wells. With this nonstationary solution of 
Eq. as we will see, the nonmonotonic evolution of the signal detection statistics of the 
bistable system can be analyzed theoretically. 



2.2 Signal detection statistics of the bistable system and the matched 
filter 

We are interested in detecting the input pulse signal s(t) from the observation of the outputs 
of bistable system or matched filter. The input pulse signal s(t) is emitted at a rate of one 
waveform every T p and lasts over a pulse duration T p . In this detection problem, we assume 
that the interval T p is both known at the emitter and the detector. But, the waveform 
amplitudes +A and —A are not known at the detector. In order to detect s(t) = +A or 

s(t) = —A, we sample the system output signals at equispaced times tj — jTp for j = 1, 2, 

The probability density p(x, T p \s(t) — ±A) of Eq. JSJ is then easily obtained for the bistable 
system. For the input bipolar pulse signals s(t) mixed with additional Gaussian white noise 
rj{t), the impulse response of matched filter is h(t) = s(T p — i), where h(t) is a scaled, 
time-reversed and shifted version of the pulse signal s(t) (0 < t < T p ). If the output of 
the matched filter is sampled at jT p that maximizes the output signal-to-noise ratio, the 
probability density p(x, T p \s(t) = ±A) has a Gaussian distribution 

p(x,T p \s(t) = ±A) = * exp[-^±^], (6) 

V27T<7 £<J 

with a = V / A 2 T P /2D [49]. 



3 



2 
1 


-1 

-2K 



(A) 



50 



150 




time t (s) 



Figure 1: Time evolution of the signals for the system of Eq. J2J) with r a = 0.1 s and Xf, = 5 
V. The sampling time step At = 0.001 s. (A) The input pulse signal s(t) is with A = 1.9 
V and T p = 1 s. s(t) = +A appears randomly in successive time intervals [(n 
(n = 1,2, . . .); (B) The mixture of signal plus noise (D = 0.1 V 2 /Hz). 



l)T p , nT p ] 



In this paper, the signal detection statistics of the bistable system and the matched filter 
are measured by the false alarm probability 

r+co 

P FA = / p{x,T p \s(t) = -A)dx, (7) 
J i 

and the detection probability 

p + ca 

Pd = / p(x,T p \s(t) = +A)dx, (8) 
J i 

where £ is the detection threshold, and the positive pulse signal s(i) = +A, shown in Fig.0 
is considered as the information carrying part of the signal in this detection problem. 

Figure (A) compares the signal detection performance of the bistable system and the 
matched filter. As the noise intensity D increases from 0.1 V 2 /Hz to 2.0 V 2 /Hz, the signal 
detection statistics of the bistable system displays a non-monotonic behavior, i.e. the SR 
phenomenon in signal detection. The detectability of the bistable system is optimized at 
D = 0.8 V 2 /Hz. The discrete points, illustrated in Fig. 0(A), represent the corresponding 
numerical simulation results of the bistable system. Here, we numerically integrate the 
stochastic differential equation of Eq. (0) using a Euler-Maruyama discretization method 
with a small sampling time step At <C r a [66]. The theoretical model of Eq. JSJ), although 
approximate, clearly describes the SR effect in the bistable system, and it agrees with the 
corresponding numerical results. Furthermore, the output signals of the bistable system and 
the matched filter are also given in Figs. El (B), (C), (D) and (E) for different noise intensities. 
The SR-type effect in the bistable system is also visible, as shown in Figs. 0(B), (C) and 
(D). 

In Fig. the bistable system, even optimized at D = 0.8 V 2 /Hz via the SR effect, can 
not outperform the matched filter for the case of Gaussian white noise. This fact has been 
pointed repeatedly by many researchers in a variety of investigations [30-32,46,50-52,61]. 
However, we have assumed complete knowledge of the probability density functions under 
hypotheses Hq and Hi, from which the conclusion of the bistable system being suboptimal 
by comparison with the matched filter is deduced. In next section, we will turn to more 
realistic problems in which some parameters of the signal are not completely known. New 
merits of the bistable system will be observed. 



4 




1 1 La^av^JLJLJ ^ 



50 100 150 



0.4 0.6 



50 100 
time t (s) 



Figure 2: (A) The detection statistics of the bistable system (solid lines) with r a = 0.1 s 
and X b = 5 V at D = 0.8 V 2 /Hz, 0.1 V 2 /Hz and 2.0 V 2 /Hz from the top down, and the 
matched filter (dashed line) at D — 0.8 V 2 /Hz. Numerical results of the bistable system are 
also presented at D = 0.8 V 2 /Hz (circles), D = 0.1 V 2 /Hz (squares) and D = 2.0 V 2 /Hz 
(down triangles), respectively. The output signals x{t) of the bistable system at (B) D = 0.1 
V 2 /Hz, (C) D = 0.8 V 2 /Hz and (D) D = 2.0 V 2 /Hz; (E) The output signal x(t) of the 
matched filter at D = 0.8 V 2 /Hz. At = 0.001 s, A = 1.9 V and T p = 1 s. 



3 Evaluating the bistable system versus the matched 
filter in unknown detection problems 

In this section, we mainly study two signal detection problems: the first is detecting the 
pulse signal with disturbed pulse rates 1/T P ; the second is a situation that the return pulse 
signal is delayed by the propagation time of the signal through the medium, resulting in the 
unknown arrival time to. Both the bistable system and the matched filter are employed as 
detectors in the two signal detection problems. The performances of the bistable system and 
the matched filter, with different measurements in different detection tasks, are compared 
thoroughly. 



3.1 The performance of the bistable system and the matched filter 
in detecting the pulse signal with disturbed pulse rates 

Consider the heuristic problem of detecting the received pulse signals with changed pulse rates 
(i.e. multi-rate pulse signals), for example, which may occur due to the desynchronization 
of devices during the detection of moving objects. For a given noisy environment, we can 
deduce the response time of the bistable system in terms of the signal amplitude A, system 
parameters r a and Xb (ref . Appendix A) . If the disturbed pulse duration T p is not too large 
(or small) with respect to the system response time T r , it will be seen in the following that 
the bistable system can detect different pulse rates from a frequency spectrum analysis. It 
is well known that intersymbol errors will appear in this situation with a matched filter [49] . 
The frequency spectrum analysis of the matched filter output also shows that too many low 
frequency components degrade the matched filter's detectability. 

In numerical simulations, the reference pulse rate is set as 1/T p = 0.2 Hz. Then, for the 
matched filter, the impulse response should be the "flipped around" version of the signal, i.e. 
h(t) = s(5—t), t e [0, 5] s. We assume that the received pulse rates of signal s(t) is distributed 
as 0.5 Hz and 0.1 Hz. Figure (A) shows the disturbed pulse signal s(t) with three rates of 
0.2 Hz, 0.5 Hz and 0.1 Hz. From the output signals of the bistable system and the matched 
filter, i.e. Figs. 0(0) and (E), the corresponding power spectra are computed and plotted in 
Figs. 13(D) and (F) respectively. The power spectrum of the bistable system output, as seen 
in Fig. 0(D), clearly illustrates three pulse rates of input signal s(t) at 1/T p — 0.1 Hz, 0.2 
Hz and 0.5 Hz. However, the frequency component at 0.5 Hz, shown in the power spectrum 
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Figure 3: (A) Input pulse signal s(t) = ±1 V with multi-rates 1/T p = 0.2 Hz, 0.5 Hz and 
0.1 Hz; (B) The mixture of input s(t) plus noise r](t) (D = 0.15 V 2 /Hz); (C) The output 
signal x(t) of the bistable system with r a = 0.1 s and Xb — 2.8 V; (D) The power spectrum 
of bistable system outputs; (E) The output signal x(t) of the matched filter. The impulse 
response of h(t) is selected in terms of the reference pulse rate of 0.2 Hz; (F) The power 
spectrum of matched filter outputs. At = 0.001 s. 



6 




50 100 150 200 250 300 350 400 



50 100 150 200 250 300 350 400 





50 100 150 200 250 300 350 400 



50 100 150 200 250 300 350 400 



Figure 4: (A) Input signal s[i] is one pulse with the amplitude 2 A — 2 V and the pulse 
duration T p = 5 s; (B) The mixture of the signal s[i] plus noise (D — 0.08 V 2 /Hz); 

(C) The output signal of the subthreshold bistable system with X\, = 2.8 V and r a = 0.1 s; 

(D) The convolution output signals computed by the matched filter, y[i] = ^ J °+ M_1 
At = 0.05 s. ' 



of the matched filter output, is weaker than the spectrum lines in the low frequency domain 
1/T P < 0.1 Hz, as seen in Fig. 03(F). This indicates that the bistable system seems to suppress 
the noise at all frequency bands, but the matched filter, more like a lowpass characteristic, 
allows too much noise at low frequencies to interfere with its frequency discrimination. 

In terms of the above comparisons, we argue that the bistable system is not so sensitive to 
the spread of the pulse rate values as the matched filter. This robust feature of the bistable 
system is not trivial in signal detection. Note that there exists an optimal pulse rate matching 
the bistable system for fixed signal amplitude and noise intensity [1,2,5], while the interest 
here is the system's robust property. 

3.2 The performance of the bistable system versus the matched 
filter in detecting the pulse signal with unknown arrival time 

In this subsection, we consider the detection problem of the received pulse signal with an 
unknown arrival time to- The hypotheses can be represented as 

H Q : x[i] = rj[i], 

H t : x\i] = s[i - i ] + r)[i], (9) 

with i = 0, 1, . . . , N — 1. Here, the observation time interval is [0, T) and N = T/At. The 
received pulse signal s[i] is with only one pulse, as shown in Fig.0](A). The pulse amplitude 
is 2 A over the interval [io, io + M — 1] (M — T p /At), but the arrival time to = ioAt is 
unknown. 

For the matched filter, the test statistic can be written as 

io + M-1 

^= ie! „,r« +1 ] ? h[i\x[i\ > 7, (10) 

i=ia 

where h[i] = s[M-l-i] for i € [0, M-l] and the decision threshold 7 < e = YhLq 1 s2 M I 49 ]- 
The arrival time of input pulses is then to = ioAt+T p (e—"f) / e. In order to find the arrival time 
to, the matched filter will compute T/At convolutions in this detection problem. Figure 0] 
(D) gives an example of the convolution output signals of the matched filter. 

Considering the symmetric characteristics of bistable systems, a bias level —A is added 
to the input x[i], i.e., x[i] = s[i — io] + n[i] — A or x[i] = n[i] — A. After the bistable system 
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Figure 5: (A) Average arrival time i e versus D; (B) Loss ratio L/K versus D. For different 
values of D, each subject is simulated with K = 80 different arrival times to- A = 1 V, 
T p — 5 s, T = 400 s and At — 0.05 s. The circle symbols represent the performance of the 
suprathreshold bistable system with Xt, = 1 V < ^/27A/2 and r a = 1 s, while the down 
triangle symbols correspond to the subthreshold bistable system with Xb = 2.8 V > \fTJAj2 
and T a = 0.1 s. The performance curves of the matched filter are plotted for the decision 
levels 7 = 200 V 2 (*) and 400 V 2 (+), with h{i) = s[99 - i], i e [0, 99] and e = 400 V 2 . 



with output y[i], illustrated in Fig. 01(C), the test statistic can be expressed as 

H(y) = y[i ] > x+, i € [0, N — M + 1], (11) 

where x+ is the positive root of equation x — x 3 /X£ + A [22]. Here, we assume that the 
bistable system needs a time interval of T r , i.e. the system response time introduced in 
Appendix A, to reach the corresponding steady-state x+. The arrival time of the pulse signal 
is then approximated as to = igAt — T r . 

There are no general theories for evaluating performance of the bistable system and the 
matched filter in this detection problem [49]. In numerical simulations, the performance of 
both detectors is quantified by the average error time 

*e = £^^, (12) 

i=i 

where K is the simulation time and j — 1, 2, . . . , K. Here, to[j] is the jth arrival time of 
the signal given by the detector. The matched filter or the bistable system, however, might 
detect a totally wrong arrival time to, which is measured by the loss ratio L/K and the loss 
time L. If we take | to[m\ — to |> T p (m = 1, 2, • • • , L) as one loss time, the average error 
time t e can be rewritten as 

' ' io\j] - f o ' 



le _ n; (13) 

3=1 

Moreover, we know the computation time T c [j] of the bistable system and matched filter in 
the jth numerical simulation. Thus, the computational efficiency of both detectors can be 
quantified by the average computation time 

f e = f;2sM. (w) 
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The plots of Figs. 0(A) and (B) show the performance of bistable systems and the matched 
filter in detecting the unknown arrival time of input signals, this is, the average error time 
t e and the loss ratio L/K. It is seen that the matched filter's performance depends on 
selecting an appropriate decision threshold 7 [49]. If the decision threshold 7 = 400 V 2 , it 
will present a higher loss ratio L/K than in bistable systems. The decision threshold 7 = 200 
V 2 makes the matched filter optimized. Two kinds of bistable systems are employed in this 
detection problem: one is subthreshold (Xb > V27A/2) and the other is suprathreshold 
(X b < V27A/2) [7,22,67]. In Fig. the average error time t e and the loss ratio L/K of 
the suprathreshold bistable system is a monotonic increasing function of the noise intensity 
D. This suprathreshold bistable system is suboptimal compared to the matched filter with 
7 = 200 V 2 , but outperforms the matched filter with 7 = 400 V 2 . For the subthreshold 
bistable system shown in Fig. 03 the average error time i e decreases as the noise intensity D 
increases, and almost approaches the performance of the matched filter (7 = 200 V 2 ) before 
the loss ratio L/K begins to increase. This non-monotonic influence of the noise on the 
performance of the subthreshold bistable system is the signature of the SR phenomenon. 

Moreover, an important merit of the bistable system is its computational efficiency. The 
average computation time of the bistable system, i.e. Eq. ijTIjl. is about T c = 0.3361 s for 
processing the data observed in the time interval [0, T] (T = 400 s). The matched filter will 
use about T c = 0.9048 s to deal with the same data. This computational characteristic of the 
bistable system maybe meaningful in the optimization of speed and efficacy of information 
processing with limited resources for data handling, storage and energy supply [29,32,53]. 

4 Conclusion 

We thoroughly evaluated the performance of the nonlinear bistable system versus the lin- 
ear matched filter in pulse signal detection problems. Based on an approximate theory of 
the probability density model, the signal detection statistics of the bistable system were 
deduced theoretically. The detection performance of the nonlinear bistable system can be 
improved by adding an optimal amount of noise, i.e. the SR phenomenon, which was also 
demonstrated by numerical experiments. The theoretical and numerical results here lead to 
the fact, emphasized by many researchers [30-32,46,50-52,61], that the bistable system is 
suboptimal to the matched filter for detecting the signal with complete known parameters. 
However, the bistable system is superior in detecting a pulse signal with unknown param- 
eters. The power spectrum analysis demonstrated that the bistable system is more robust 
than the matched filter for discriminating the disturbed pulse rate. In the case of detecting 
the unknown arrival time of the pulse signal, the bistable system is also comparable to the 
matched filter. Furthermore, the detection scheme of the bistable system is more efficient 
with computational saving. These significant results verified the potential applicability of 
nonlinear bistable systems in the signal detection field, which might deserve further studies, 
for example, in exploiting other nonlinear systems exhibiting SR effects in practical detection 
problems. 
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A System response time and nonstationary probability 
density model 

A.l System response time 

In Eq. (0), the Fokker-Planck operator Lpp = -§^V (x) + ^--§^r is not a Hermitian operator 
[65]. We rescale the variables as 



X b = A = A/^Dj7 a , t = t/r a , y = x/^Dj7 ai (15) 

Eq. J2J becomes 

where V'(y) = —y + y 3 /X b 2 =p A. The steady-state solution of Eq. (jlf)fl is given by 

P {y) = lim p(y,r) = Cexp[-V(y)], (17) 

r— >oo 

where C is the normalization constant. A separation ansatz for p(y,r) [65], 

V(v) 

p(y, t) = u{y) exp[ p] exp(-Ar), (18) 

leads to 

Lu = -Aw, (19) 

with a Hermitian operator L = (d 2 /dy 2 ) — [\V 2 {y) — \V"{y)\- The functions u(y) are 
eigenfunctions of the operator L with the eigenvalues A. Multiplying both sides of Eq. I|19|) 
by u{y) and integrating it, yields 

A = J-™W 2 (V) + u 2 (y)[\V 2 (y) - \V"{y)]}dy 

where eigenfunctions u(y) satisfy the boundary conditions of \im y —,± 00 u(y) — and lmiy^± 00 u' 
0. The eigenvalue problem of Eq. H19[) is then equivalent to the variational problem consisting 
in finding the extremal values of the right side of Eq. (120ft . The minimum of this expression is 
then the lowest eigenvalue Ao = 0, corresponding to the steady state solution of Eq. I|17fl [65]. 
Here, we adopt eigenfunctions u(y) — p(y) exp[— V(y)/2] andp(y) ^ 0, Eq. (|2*U)l then becomes 

x _ f-^{P /2 (y) + \p 2 {y)V l2 (y) \[V'{y)p 2 {y)]'}ex V [~V{y)]dy 
x ~ p+oo T, ; r 77 , • \ ll > 



J-™P 2 (y) exp[-V(y)]dy 



Since 



[V'(y)p 2 (y)]'exp[-V(y)}dy = V (y)p 2 (y) exp[-V (y)}\ 



+ OC 



r+oo 

p 2 (y)V 2 (y)exp[-V(y)]dy = / p 2 (y)V' 2 (y) exp[-V(y)]dy, 



Eq. l(2*T)) can be rewritten as 

x f+™ P > 2 (y)eM-v(y)}dy 

f-™P 2 (y) exp[-V(y)]dy 
Assume p(y) = do + d\y + ■ ■ ■ + d n y n and the order n is an integer, we obtain 

([K]-\[M]){d}=Q, (23) 
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with eigenvectors {d 1 } = [d l , d\,. . . ,d l n ] corresponding to eigenvalues Aj for i = 0, 1, . . . , n. 
The elements of matrices [M] and [K] are 

/+oo 
y i+j exp[-V(y)]dy > 0, 
-OO 

/+oo 
ijy l+J - 2 eM-V(y)}dy>0, 
-oc 

where i, j = 0, 1, . . . , n. The matrix [M] is positive definite and the matrix [K] is semi-positive 
definite. The minimal eigenvalue Ao is zero. The inverse of minimal positive eigenvalue Ai 
describes the main time of the system tending to the steady state solution of Eq. p7fl . what 
we call the system response time. Note the time scale transformation in Eq. (JT5J, the minimal 
positive eigenvalue should be Ai/r a and the real system response time is T r = T a j\\. 



A. 2 Nonstationary probability density model 

From Eq. Ij23(l . we can obtain the eigenfunctions Ui(y) = Pi(y) exp[— V(y)/2] corresponding 
to the eigenvalue Ai for i = 0, 1, . . . , n, where Pi(y) = d l Q + d\y + • • • + d l n y n . The eigenvectors 
{d 1 } = [d l , d\, . . . , d l n ] are normalized. Because L is a Hermitian operator, eigenfunctions 
Ui(y) and Uj(y) are orthogonal 

/+oo 
Ui(y)uj(y)dy = 6ij, (24) 
-OO 

where i,j = 0,1,..., n. p{y,r) can be expanded, according to eigenfunctions Ui{y) and 
eigenvalues A 2 ;, as 

(0 (y, T ) = ^C iUi ( 2 /)e X p[-^M]e X p[-A i r], (25) 

i=0 

where Ci are normalization constants. In this paper, we take an approximate expression of 
p(y,r) ~ Yli—o Ci u i(y) ex P[ — V(y)/2] exp[— Aj-r] instead of Eq. J2HJ). Then, if the preceding 
input signal is s(t) = ^fA and the next one is s(t) = ±A, a simple nonstationary probability 
density model is derived as 

p(x, t\s(t) = ±A) ~ p(x\s(t) = ±A) + [p(x\s(t) = tA) - p{x\s(t) = ±A)} exp(-t/T r ), (26) 

with the initial condition p(x,t = 0|s(t) = ±j4) = p(x\s(t) — ^A) and the stationary condi- 
tion p(x,t = +oo\s(t) = ±A) = p(x\s(t) = ±A), respectively. Here, p(x\s(t) = ±A) are the 
steady-state solutions given in Eq. QJ, and V(x) = — x 2 /2 + x 4 / (4X£) =p Ax correspond to the 
constant inputs s(t) = ±A respectively. This kind of nonstationary solution of p(x,t\s(t) = 
±A) can be further developed into the case of p(y 7 r) ~ J27=o CiUi(y) exp[— V(y)/2] exp[— A^r] 
for n > 2. 
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